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Abstract 

It has been shown that a final state interaction model, used to describe 
J ftp suppression, can also describe the fixed pt suppression of the tt° (and 
charged pions) yield at all values of pt, with a final state interaction cross- 
section a close to one milibarn. We propose an extension of the model to the 
pion motion in the transverse plane - which introduces a dependence of the 
suppression on the azimuthal angle 9r. Using the same value of a, we obtain 
values of the elliptic flow i>2 close to the experimental ones, for all values of 
Pt, including the soft pt region. 
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1 Introduction 



Statistical QCD (i.e. QCD applied to a system in thermal equilibrium) on a 
lattice predicts the existence of a new phase of matter, quark matter or quark gluon 
plasma. It is therefore most important to find signals of the thermalization in the 
heavy ion data. To achieve thermalization is by no means trivial. High densities are 
only a necessary condition. Strong final state interaction is also needed, as well as 
an interaction time long enough for the system to equilibrate. 

Nuclear collisions are, a priori, not favorable to produce a thermalized system. 
Indeed, pA collisions can be described with independent string models [1]. There 
is no indication in the data of final state interaction between particles produced in 
different strings - despite the fact that the proton collides with several nucleons of 
the nucleus and densities several times larger than the proton proton one are reached 
and that several strings are produced in a transverse area of 1 fm 2 . 2 In a central 
AA collision this string density is higher and "cross-talk" between different strings 
is expected. Correspondingly, final state interaction is seen in the data : strangeness 
enhancement, J/ip suppression, large px suppression, elliptic flow, etc. cannot be 
described with independent strings. It could happen that string interaction is so 
strong that the very concept of string becomes meaningless. On the contrary, in 
view of the nuclear transparency observed in pA, one could expect that final state 
interaction will take place with a comparatively small cross-section. In this case the 
bulk of the produced particles is hardly affected by the final state interaction. Only 
rare events, of the type mentioned above, are strongly modified. 

Heavy ion data lend support to this scenario. For instance, one could expect that 
thermalization strongly affects the distribution of produced particles - especially 
at low pt- Experimentally, the variation of < pr > between the most peripheral 
and the most central collisions is quite small (about 20 % at RHIC) and consistent 
with the small increase observed in pA. Also the RHIC data on pjhr are flat 
from peripheral to central collisions - showing no sign of strong pp annihilation 
which would produce a change in the pjhr ratio for (out of equilibrium) peripheral 
collisions followed by a flatness when thermalization is reached. Turning to rearer 
2 This is probably due to the small transverse size of the string as characterized by the correlation 
length of the non-perturbative interaction (about 0.2 fm from lattice data). 
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events, the centrality dependence of the ratios B/h~ and B/h~ at RHIC shows no 
sign of saturation - increasing monotonically from peripheral to central collisions 
for A, S and Q - contrary to statistical model predictions [2]. Also the decrease 
of < p\ > of J/ip at large centralities, predicted in a deconfining phase transition 
scenario, is not seen in the data [3]. 

One of the main arguments for (early, transverse) thermalization at RHIC is the 
description of elliptic flow data with an ideal (no viscosity) hydrodynamical model. 
However, this view has recently been challenged [4]. Experimentally [5], it turns 
out that ideal hydrodynamics describes the v 2 data only in a small low-p? interval 
and a very small rapidity range. Moreover, it only describes the minimum bias data 
and fails to reproduce the centrality dependence. In particular, it overestimates the 
data at small impact parameter, where the densities are higher. A review of hydro- 
dynamic model [6], coalescence model [7] and transport model [8] calculations and 
its comparison with experimental data, as well as a comprehensive list of references 
on these models can be found in [5]. 

Quantitative support to the idea of a weak final state interaction, which affects 
only rare events can be found in refs. [9] and [10]. In [9] it has been shown that 
the strangeness enhancement data [11] can be described with a final state interac- 
tion cross-section a = 0.2 mb. In [10], the J/ip suppression at CERN-SPS [12] is 
reproduced with a = 0.65 mb. More recently, it has been shown [13] that the large 
p T suppression of 7r°'s, observed at RHIC can be reproduced with the same final 
state interaction model [9] [10], using an effective cross-section of the order of 1 mb. 
In this formalism the ratio Raa{Pi Pt) = dN AA / dy(b, p T ) / n(b)dN pp / dy (where n(b) 
is the average number of binary collisions at fixed impact parameter b) can be re- 
produced at all values of - including the soft pt region. Moreover, the large 
Pt results turn out to be rather insensitive to the value of the pt shift of the 7r° 
resulting from its interaction with the hot medium. (In our approach the results for 
charged pions and 7r°'s are identical). 

In this paper we propose an extension of the final state interaction model [9, 10, 
13] in order to take into account the motion of the pion in the transverse plane. This 
introduces a dependence of the pion survival probability on its azimuthal angle 9 R , 
which, in turn, gives rise to a (positive) contribution to the elliptic flow i> 2 . 
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The fact that the mechanism responsible for the large pr suppression does give 
a contribution to the elliptic flow at large pt of the order of the experimental one is 
known [14] [15]. Here we show that, in the model [9, 10, 13], values of v 2 very close 
to the experimental ones are obtained in the whole pt region, including the low pr 
one, using the same value of the final state interaction cross-section (of the order of 

1 mb) needed to reproduce the experimental values of Raa- 

Although this contribution to i>2 results from an asymmetry in the azimuthal 
angle, it can be qualified as non-flow [15]. Indeed, the mechanism from which it 
arises (fixed pt suppression) is maximal at zero impact parameter and, moreover, 
thermalization is not needed. 

2 The final state interaction model 

As in refs. [9, 10, 13] the model we are going to use is based on the gain and 
loss differential equations which govern final state interactions. Following [16] we 
assume boost invariance and dilution in time r of the (transverse) densities pi due 
only to longitudinal expansion (p ~ 1/r). We then have [16] 

r = Yl °m Pk Pt ~ a ik Pi Pk ■ (1) 

Here pi = dN AA ~* l (b) / dyd 2 s are transverse densities and are the final state inter- 
action cross-sections averaged over the relative velocities of the incoming particles. 
The first term of (1) describes the gain in type % particle yield resulting from the 
interaction of k and t. The second one corresponds to the loss of type i particles 
resulting from its interaction with particle k. Eqs. (1) have to be integrated between 
initial time r and freeze-out time Tf. The solution depends only on the ratio t//t . 
Actually particles i,j,k,£--- can be either hadrons or partons. Indeed, at early 
times the densities are very high and hadrons not yet formed. Thus, at early times 
eqs. (1) describe final state interactions at a partonic level. Only at later times, we 
have interactions of full fledged hadrons and, thus, a represents an effective cross- 
section averaged over the interaction time. For this reason formation times - which 
are not introduced here - do not play an important role. 

Let us now consider a n° produced at fixed px interacting with the hot medium 
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consisting of charged plus neutrals. We assume that in the interaction, with cross- 
section a, the 7i° suffers a decrease in its transverse momentum with a pr-shift 5pr- 
This produces a loss in the n° yield in a given p T bin. There is also a gain resulting 
from 7r°'s produced at pr + Spx- Due to the steep fall-off of the px spectrum the 
loss is larger than the gain, resulting in a net suppression of the 7r° yield at a given 
p T . In this case there is a single final state interaction channel and eqs. (I) can be 
solved analytically. Integrating both members of eq. (I) over d 2 s, we obtain the 
following expression for the survivable probability of the tt° in the medium 

N n o(p T + 5pr) 



Sno (y, pr, b) = exp < -a p(b, y) 



£n(p(b,y)/p pp (y))\ ■ (2) 



nApt) 

Here N n o(p T ) denotes the yield of 7r°'s at a given p T (see section 3). With 5p T — > oo, 
the quantity inside the square bracket is equal to one and the gain term vanishes. 
In this case, the survival probability has the same expression as in the case of J ftp 
suppression without cc recombination. With Spx = 0, the loss and gain terms are 
identical and the survival probability is equal to one. 
The density p(b, y) in eq. (2) is given by 

dN AA 

P (b,y) = —(b)/G(b) (3) 

where the inclusive AA distribution refers to the production of charged plus neutrals 
and G(b) is an effective transverse area (see below). 

The argument of the logarithm in eq. (2) is the ratio of the freeze-out 77 to the 
initial time r . The solution of (I) depends only on the ratio T//r . We use the 
inverse proportionality between time and density (see above) to write 

r/ /r = po/pf = p(b,y)/p PP (y) (4) 
where the freeze-out density pf has been taken to be the pp one, with 

P P M = 3 I nRl . (5) 

At y/s = 200 GeV, p pp (y = 0) = 2.24 fm~ 2 . 3 

3 With this prescription the average value of 17/70 is of 5 -f- 6. Thus, if to ~ 1 fm we obtain a 
freeze out time 77 ~ 5-6 fm. 
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The densities in eqs. (l)-(4) are computed in the dual parton model [1]. They are 
given by a linear combination of the average number of participants and of binary 
collisions. The coefficients are obtained from convolutions of momentum distribution 
and fragmentation functions. Their numerical values for Au Au collisions at y/s = 
200 GeV at each b are given in [17]. Shadowing corrections, described in detail in 
[17], are included in the calculation. Finally the geometrical factor G(b), resulting 
from the integration over d 2 s of the two members of eqs. (1), is given by 

r dN AA I dN AA 

where n(b, s) = a pp TA(s)TA(b — s)/aAA{b). The nuclear profiles Ta are computed 
from Woods-Saxon densities [18]. We use the non-diffractive inelastic cross-section 
<jpp = 34 mb at y/s = 200 GeV. In writing eq. (6) we have assumed that the large 
Pt probe scales with the number n of binary collisions 4 . 



3 Numerical results for the tt° suppression 

The suppression of the 7r° yield at fixed px is usually presented via the ratio 

Raa(P,Pt) = —^-(P'Pt) / < h ) -fy-M ■ ( 7 ) 

In this paper we shall concentrate on the mid-rapidity region \y*\ < 0.35. In order 
to compute this ratio we have to know its value R AA (b,pr) in the absence of final 
state interaction. Its pt integrated value (close to 1/3 for the 10 % most central 
Au Au collisions) is obtained in the dual parton model with shadowing corrections 
[13]. Its pt dependence (p? broadening) is controlled by two mechanisms. The 
first one is the decrease of shadowing effects when p T increases. This produces an 
increase of R AA (b,pT) towards unity. The second mechanism is the Cronin effect 
proper, which produces a rise of pr) above unity. In ref. [13] we have used a 

phenomenological parametrization of both effects of the form 



4 In the dual parton model at very high energies, dN AA /dyd 2 s, before shadowing correc- 
tions are introduced, is also proportional to the number of binary collisions. In this case, 
C7 _1 (fe) = J d 2 s{TAA(b, s)] 2 /T AA (b) is a purely geometrical factor. The results obtained with these 
two definitions of G(b) are quite similar. 
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where po(b) = (n — 3)/2 < px >b, n = 9.99 and < px >b is the experimental 
value of < px > at each b [19]. Details of the phenomenological parametrization of 
the pt distribution leading to eq. (8) can be found in [13]. The initial suppression 
R AA (b,px) obtained from (8) reaches values substantially larger than unity and does 
not decrease towards one at large px (see Fig. 1 ; numerical values are given in [13]). 

The value of the n° suppression, characterized by the ratio (7), can be com- 
puted by multiplying R AA (b,px) by the n° survival probability due to final state 
interaction, given by eq. (2): R AA (b,p T ) = R° AA (b,p T ) S n o(pr,b). 

In order to compute the survival probability, we have to define the quantity 
inside the square brackets in eq. (2). Using the same parametrization of the pt 
distribution as in eq. (8) we have 

N n o(p T + 5p T ) = f pojb) +p T + Sp T \~ n /q\ 
NApt) U I Po(b)+p T J • U 

Note that this quantity represents the gain of particles at a given p? resulting 
from the particles or partons produced originally at px + 5px and that are shifted 
due to the interaction with the dense medium. As for eq. (8), this parametrization 
is extracted from the power-law function proposed in Ref. [19]. 

We have studied different possibilities for the shift 5p T [13]. In a way, this is 
equivalent to a test on how the mean energy loss (px being E at y = 0) behaves 
with the energy. We have also tried different parametrizations for R AA (b,px), i.e. 
the ratio in the absence of final state interaction. We have found that our final result 
depends little on the form of R AA when we take a p^-shift of the form Spx = p T /C. 
We have also check that our results for px > 5 GeV/c are almost independent on 
the form of the shift. 

A reasonable description of the data [20] for all values of pr has been obtained in 
[13] using 

Sp T = p% 2 J (20 GeV 1/2 ) (10) 

and the initial suppression R AA given by eq. (8). Note that, at low p?, 5pr has to 
decrease in order to match the dual parton model results. For obvious reasons Spx 
is expected to vanish at small pt- 

This pt shift can be used in the whole range of pr available. However, as men- 
tioned above, one can also use a constant Spx for p T > 6 4- 7 GeV, without any 
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significant change in the results [13]. Using (10) in the whole pt region and a value 
of a = 1.3 mb, we obtain the results shown in Figs. 1 and 2. In Fig. 1 we show 
our results for R AA (p T ,b), eq. (7), in the centrality bin 0-10 %, using a = 1.3 mb. 
They agree quite well with the available data [20] in the whole pr region, including 
the soft one. The corresponding values of the initial suppression R AA (p T , b) are also 
shown. In Fig. 2 we present our results for the dependence of Raa on the number 
of participants, N part , for pt > 4 GeV. The comparison with the data shows that 
our model reproduces the centrality dependence. 



4 Extension of the model to the transverse 
motion 

In its formulation, the final state interaction model introduced in section 2 takes 
into account the longitudinal expansion - with no consideration for the motion in 
the transverse plane. Elliptic flow, on the contrary, results from an asymmetry in the 
azimuthal angle, and, thus, the motion in the transverse plane plays a fundamental 
role. The extension of the model to take this transverse motion into account is by 
no means trivial. To the best of our knowledge a totally satisfactory formulation is 
not available. In the following, we propose a simple extension of the model taking 
only into account the different path length of the 7r° in the transverse plane for each 
value of its azimuthal angle Or - measured with respect to the reaction plane. At 
y* ~ 0, the path length Re R , measured from the center of the interaction region 
(overlap of the colliding nuclei) is given by 

sin V R 

where R A = 1-05 A 1 / 3 fm is the nuclear radius and sin a = bsm.9n/2R A . Note that 
eq. (11) is only valid in the region < 9r < 90°, and the integral in 6r from to 
360° is obtained by integrating from to 90° and multiplying the result by four. 
Our ansatz consists in the following replacement in eq. (2) 

p(b, y*~0)^ p(b, y* ~ 0)R e J < R 0R > (12) 

where < Rq r >= J 90 d0 R R Sr J J 90 d0 R . This is motivated by the fact that, the 
duration of the interaction, characterized by the argument of the logarithm in eq. 
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(2), as well as the density of the medium traversed by the 7r°, are expected to be 
proportional to the n° path length associated to its azimuthal angle 9r, inside the 
overlap region of the colliding nuclei. Due to the division by < Rg R > in (12), the 
results in Section 3 are practically unchanged. 

With the replacement (12), the survival probability of the 7r°, S^ipr, b), depends 
on the angle 9 R and the elliptic flow can be obtained as 

, i> C'de R cos20 R SlSiPT.b) ,,,, 

M= c^s: siPT , b) ■ (13) 

Clearly, when the 7r° moves along the reaction plane its path length will have its 
minimal value and the survival probability its maximal one. On the contrary, for 
Or = 90° the path length will be maximal and the survival probability minimal. 
Thus, the resulting value of v 2 is positive. 

Using the same value a = 1.3 mb of the final state interaction cross-section and 
the same pr-shift introduced in section 3 in order to describe the experimental values 
of RAA(b,pr), we obtain the values of v 2 versus pr and centrality shown in Figs. 3, 

4 and 5. These results are compared with PHENIX and STAR data [21, 22, 23] 
for charged hadrons. Notice that the results for identified pions are practically the 
same at least in the low pr range [5]. 

It is quite satisfactory that our results are in agreement with experiment for 
central and medium central collisions, which contain the bulk of the data. The 
situation for peripheral collisions is further discussed in Section 5. We also describe 
v 2 in the low and intermediate pt region. At pr > 3 GeV our results (dashed line) 
overestimate the data. A natural way to overcome this problem is discussed in the 
next section. 

5 The large £>t region 

Using eq. (8), we have obtained an initial suppression Raa(Pt) that increases 
rapidly towards unit with increasing p^, reaches values substantially larger than 
unity and does not decrease towards one at large pt (see Fig. 1). Such a behaviour 
is expected in the small pr region where R AA increases rapidly. As mentioned above 
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part of this increase is due to the decrease of shadowing when increases and part 
of it to the Cronin effect. At large pr, on the contrary, this behaviour is not the 
conventional one. One rather expects that at large p T , say px > 5 GeV, where the 
shadowing corrections are already quite small, the value of R AA is close or equal to 
unity - consistent with the smallness of the Cronin effect observed in d-Au as well 
as in the direct photon production data in Au Au. 

The value of the n° suppression, characterized by the ratio (7) is computed by 
multiplying R AA (b,pr) by the n° survival probability S n o due to final state inter- 
action, given by eq. (2). In view of that, it is not possible to determine R° AA from 
the large px suppression data alone - since one can change the two factors R AA 
and S n o, keeping their product unchanged. On the contrary, we have seen that the 
values of the elliptic flow t> 2 depend only on S n o. Thus, the combined analysis of 
both sets of data allows a separate determination of the two factors. We are going 
to show that only values of R° AA close to unity at large p T can describe both sets of 
data. Accordingly, we proceed as follows: we use eq. (8) only for small pr and take 
R AA (b,pT > 5 GeV) = 1. Since the large px suppression is given by the product 
Raa(Pt) S-k°{Pt) , this change of R AA for p T > 5 GeV has to go along with a change 
in S n o(p T ) - which can be achieved by changing the pr-shift at large pr- Let us 
parametrize it in the form 5pr = Pt/C with the parameters a and C determined 
in such a way that the size and (practically flat) shape of the large px suppression 
data is preserved. This is achieved with a = 0.8 and C = 9.5 GeV -0 - 2 . Thus, we 
shall take 

S PT = p° T 8 /9.5 GeV- - 2 for p T > 5 GeV . (14) 

Let us note that eqs. (10) and (14) lead to the same value of the p^-shift at 
Pt = 2.9 GeV. Therefore, in order to have a smooth p^-shift in the whole px region, 
we are going to use eq. (10) for p T < 2.9 GeV and eq. (14) for p T > 2.9 GeV. In the 
first region (p T < 2.9 GeV) we use the form of the initial suppression R AA (p T ) given 
by eq. (8) and, at large p T) we take R AA (p T > 5 GeV) = I. In this way our results 
for R AA are practically the same, as can been seen in Fig. 1, where we show the form 
of the two different parametrizations for the ratio R AA - ratio before the interaction 
- and the resulting Raa- Both parametrizations give similar results concerning the 
ratio for large px, for this reason Fig. 2 remains practically unchanged. However, the 
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values of v 2 at large px are substantially reduced. One obtains the full line in Fig. 3. 
This shows a saturation of V2 at pr > 3 GeV which is consistent with experiment. 

More information on the large pt values of v 2 and their centrality dependence is 
provided by the preliminary data on the azimuthal angle dependence Sb r (pt) of the 
7T° suppression at different centralities [24] . Our results are given in Fig. 6 and com- 
pared with experiment. The agreement with the data is reasonable, indicating that 
the centrality dependence is approximately reproduced. This is in sharp contrast 
with the disagreement observed at lower values of pr (Fig. 4), where our centrality 
dependence is too weak as compared with the data but in much better agreement 
with the ones in Fig. 5. This important point has to be clarified. 

In the above calculation of V2 we have evaluated the path length of the 7r° (eq. 
(11) at its most probable position in the transverse plane, namely the center of 
the interaction region determined by the overlap of the two colliding nuclei. In a 
more accurate calculation, we should evaluate this path length at the point in the 
transverse plane where the 7r° has been produced. However, we have estimated 
that the error induced by this approximation is less than 10 %. Indeed, we have 
checked numerically that moving the 7r° production point along the Or = direction 
the value of i>2 decreases, while, moving it in the Or = 90° one, it increases. The 
average of these two values of i>2 is close to the one obtained using (11) - although 
somewhat smaller. The difference increases when one moves away from the center 
and reaches the 10 % level at the mid points on the two axis. Due to the integration 
over TAA(b,s) present in our calculation, the average position of the produced n° 
is closer to the center of the interaction region than to the mid-points. Thus, the 
estimated error is less than 10 %. 

While the approximation discussed above overestimates slightly the value of V2, 
there is a physical effect which can increase its value in a more substantial way. 
Namely, if, besides the pr-shift, the final state interaction produces also a shift 50r 
in the azimuthal angle Or of the 7r°, it can be shown that the value of v 2 increases, 
irrespective of the sign of 50 R . For instance with \S0 R \ = 15°, the increase of v 2 at 
b = 6 fm and pr = 1.35 GeV can be as large as 30 %. 5 

5 The effect depends on the time at which the shift 59r takes place. It is maximal when this 
time is close to freeze-out. 
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6 Conclusions 



We have proposed a final state interaction model which takes into account 
the different path length of a particle in the transverse plane for each value of its 
azimuthal angle in the overlap area of the colliding nuclei. The model is an extension 
of one previously introduced to describe stangeness enhancement, J/ip suppression 
and fixed px suppression. Using the values of the final state interaction cross- 
section (about 1 mb) and of the pr-shift which describe the RHIC data on fixed pt 
7t° suppression, we obtain values of the elliptic flow v 2 close to the measured ones 
for all values of p T) including the soft region. As for the centrality dependence, the 
values of v 2 are well reproduced in the region b < 6 fm which contains the bulk of 
the data. We argue in Section 5 that the situation for peripheral collisions has to 
be clarified. 

Moreover, combining the analysis of large pt suppression and elliptic flow data, 
it is possible to obtain information on the initial suppression, R aa (pt), i.e. the 
one in the absence of final state interaction - which results from the combined 
effect of shadowing and the Cronin effect. Values of R AA close to unity at large p T 
(Pt <L 5 GeV) are required in order to describe both data sets. 

The mass dependence of v 2 is reproduced quite well. The decrease of v 2 with 
increasing mass is just a consequence of the broadening of the pt distributions when 
the mass increases. It will be described in detail in a forthcoming publication [25]. 

In view of the uncertainties discussed at the end of section 5 we have restricted 
our calculation to the second order coefficient v 2 and to mid-rapidities. However, 
our formalism allows to understand qualitatively the strong decrease of v 2 when 
moving away from mid-rapidities. This decrease follows the decrease of the medium 
density as expected in our approach. Thus, we expect a dramatic decrease of v 2 
when moving away from mid-rapidities. 

We do not claim that the mechanism we have introduced gives the only contribu- 
tion to the elliptic flow. The latter is a very subtle observable to which presumably 
different mechanisms do contribute. In our opinion, our knowledge of the dynamics 
of the nuclear interaction is not sufficient to disentangle all these mechanisms and, 
therefore, to allow to draw clear-cut conclusions regarding the interpretation of the 
measured values of v 2 . 
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Together with previous work on strangeness enhancement, J/ip suppression and 
fixed pt suppression, all of which can be described in our framework with final 
state interaction cross-section smaller than or of the order of 1 mb, this work lends 
support to the idea that, despite the large densities reached in central Au Au colli- 
sions, the final state interaction is rather weak. In terms of string models, it means 
that there is "cross-talk" between different strings. However, the concept of string 
remains useful and these models can be used to compute the densities needed as 
initial conditions in the gain and loss differential equations which govern the final 
state interaction. Such a weak cross-talk, in the presence of many strings per unit of 
transverse are, is supported theoretically by the small transverse size of the string - 
with a radius of the order of 0.1 fm. On the experimental side it is supported by the 
nuclear transparency of proton-nucleus collisions - where no "cross-talk" is needed 
to reproduce the data. As a consequence, only rear events are substantially affected 
but the bulk of the system is not - indicating that the strength and duration time 
of the final state interaction are not enough to drive the system to thermal equili- 
bration. 
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Figure 1: The tt° suppression factor RauAu(pt, b) at ^/s = 200 GeV, eq. (7), versus 
Pt i n the centrality bin 0-10 % and the corresponding tt° initial suppression fac- 
tor RauAu(Pt, b) (i.e. in the absence of final state interaction). The dashed lines 
correponds to our results obtained using eq. (8) for R\ u Au an d a pr-shift given by 
eq. (10). The continous lines are obtain with R° AuAu described in Section 5 and the 
Pr-shift given eq. (10) for p T < 2.9 GeV and eq. (14) for p T > 2.9 GeV. The data 
are from [20]. 
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Figure 2: The n° suppression factor RauAu(pt, b) versus the number of participants, 
Npart, for pt > 4 GeV. The data are from [20]. 
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Figure 3: Values of the elliptic flow v 2 (b,p T ) versus p T in the centrality bin 13 %- 
26 %. The continous line correponds to our results obtained using eq. (14), the 
dashed line corresponds to our results obtained using eq. (10). The data are from 
[21]. 
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Figure 4: Values of the V2(b, pr) versus the number of participants for different values 
of pt ■ Pt = 0.4 GeV (lower line), px = 0.75 Gev (middle) andpr = 1-35 GeV (top). 
The black data are from PHENIX [21], the open data are from STAR [23]. 
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Figure 5: Values of the v 2 (b, px) versus the number of participants for different values 
of pt ■ Pt = 0.4 GeV (lower line), px = 0.75 Gev (middle) and pr = 1.35 GeV (top) 
compared to: data in minimum bias collisions from PHENIX [22] (black symbols) 
and STAR [23] (opened symbols) -above-, data from PHENIX [22] -belov^. 
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Figure 6: Values of the n° suppression RauAu{P^6r) as a function of the azimuthal 
angle 9r (measured from the impact parameter direction) for p T > 4 GeV in various 
centrality bins. The continous line correponds to our results obtained using eq. (14), 
the dashed line corresponds to our results obtained using eq. (10). The preliminary 
data are from [24]. 
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